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ABSTRACT 

We model a cold-wall atmospheric pressure impinging jet reactor 
to study the CVD and gas-phase nucleation of Ti0 2 from a titanium 
tetra-iso-propoxide (TTIP) /oxygen dilute source gas mixture in 
nitrogen. The mathematical model uses the computational code FIDAP 
and complements our recent asymptotic theory for high activation 
energy gas-phase reactions in thin chemically reacting sublayers. 
The numerical predictions highlight deviations from ideality in 
various regions inside the experimental reactor. Model predictions 
of deposition rates and the onset of gas-phase nucleation compare 
favorably with experiments. Although variable property effects on 
deposition rates are not significant (-11% at 1000K) , the reduction 
of rates due to Soret transport is substantial (-75% at 1000K) . 


1. INTRODUCTION 

Production of sophisticated materials with superior properties 
by CVD requires an understanding of coupled transport and chemical 
processes. This goal can be realized in research programs by 
combining both experiments and modeling. 

Recent Yale University research focused on TiO a CVD with simul- 
taneous gas-phase reaction leading to particle nucleation [1-2]. 
A simple asymptotic theory treated the onset of reactions leading 
to particle nucleation and reduced CVD rates at high surface tempe- 
ratures. The theory assumed that all gas-phase reactions are rest- 
ricted to a thin chemical sublayer adjacent to a hot CVD surface 
due to their large activation energies. They also developed a lab- 
scale, cold-wall, impinging jet CVD reactor [3]. Using experiments 
and appropriate asymptotic and numerical models, we expect to guide 
future reactor designs and scale-up, resulting in maximum CVD rates 
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while avoiding harmful particle nucleation. 

The present axisymmetric, numerical model points to the role of 
bouyancy and recirculation in the experimental reactor and clari- 
fies the effect of Soret transport in modifying TTIP mass transfer 
rates. It provides the ability to capture the transition from 
surface kinetics to gas-phase transport control. Furthermore, it 
extends the earlier asymptotic theory [1] by removing the restric- 
tion that the activation energies of gas-phase reactions should be 
"large" (i.e. the chemical sublayer should be "thin"). 

For the case studied at Yale [1-3] for the CVD of Ti0 2 from di- 
lute TTIP/0 2 in N 2 , the ratio of chemical sublayer to thermal bound- 
ary layer thickness is estimated to be 1/5 at the onset of particle 
nucleation (T=1050K) , which may be too large for the accuracy of 
the asymptotic theory. Hence, our numerical model may ultimately 
be better able to capture the effects of homogeneous chemistry 
leading to an observed deposition rate fall-off above -1050K. 


2 . EXPERIMENT 

We use a cold-wall, atmospheric pressure, axisymmetric, imping- 
ing- jet reactor (Fig. 1) . The liquid TTIP source is a constant 
temperature bubbler. After mixing TTIP vapor with excess oxygen 
(to "burn away" any co— deposited carbon) and diluting with nitrogen 
in a short mixing chamber, the gas jet emerges from a converging 
nozzle and impinges on a polished quartz substrate (diam. 1.3cm). 
The substrate and the alumina substrate holder are supported from 
below by an RF-heated graphite susceptor. The susceptor tempera- 
ture is measured using a Pt-30%Rh/Pt-6%Rh thermocouple just below 
the substrate. This thermocouple reading is correlated with direct 
substrate surface temperature measurements made in situ using known 
melting point lacquers (with -3-4% accuracy) . Gas and surface tem- 
peratures at several other locations are also measured by thermo- 
couples. The concentration of water vapor in the reactor is less 
than 3ppm. Deposition rates are measured by in situ interfero- 
metry and confirmed by ex situ weight gain. Further details of the 
experimental system and operating procedure are given in [3]. 

3. PHYSICO-CHEMICAL/NUMERICAL MODEL 

Our modeling study adopts the finite-element-based computational 
fluid dynamics code FIDAP. The new version of FIDAP incorporates 
many phenomena relevant to CVD processes, such as temperature 
dependent fluid density, transport and thermodynamic properties, 
Soret diffusion, and gas-phase and surface chemical reactions. 

We approximate the reactor by using a 2-D axisymmetric geometry 
(Fig. 2) . The inlet and exit are placed sufficiently far from the 
deposition zone to eliminate numerical uncertainties on rate 
predictions. The numerical mesh above the substrate is fine enough 
to resolve the chemical sublayer. Temperature boundary conditions 
are interpolations of measured surface temperatures in the Yale 
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axi symmetric impinging jet reactor. 

We assume that the gas-phase decomposition of TTIP is first or- 
der with a rate constant of 3 . 96xl0 5 exp (-8480/T) s' 1 (T is in Kel- 
vins) [4]. We also assume that the gas-phase products form nonde- 
positing particles. We fit our low substrate temperature rate data 
to a Soret-corrected Arrhenius expression, 1 . 2lxl0 9 exp (-16480/T) 
m/s, to obtain a pseudo-first order surface reaction rate constant. 
This expression corresponds to a TTIP reactive sticking probability 
of unity at 987K. For simplicity, the surface reaction rate 
constant is kept at its 987K value for higher substrate tempera- 
tures so that the sticking probability does not exceed unity. 

The Soret diffusion factor for dilute TTIP in N 2 is calculated 
from kinetic theory (the estimated Lennard-Jones parameters for 
TTIP are a=8.13A and e/k=589K) and fitted to the expression 
a T =l . 97 i [ i- (223 . 6) /T] . Unfortunately, FIDAP is currently 
restricted to constant Soret dif fusivities , Dj = pDci)Q: T , hence we 
evaluate the gas (N 2 ) density p, TTIP Fickian diffusivity D, and a T 
at T fl]m = (T nozzle +T surface ) / 2 , and use the TTIP mass fraction a at the T fllm 
location for the corresponding diffusion limit calculation. 

4. RESULTS AND DISCUSSION 

First, we compare FIDAP predictions to the observed flow field 
inside the reactor. The experiments use argon seeded with fine ti- 
tania particles for flow visualization. Figure 3 shows a typical 
comparison at a Reynolds number of 750 based on nozzle diameter at 
STP and 900K substrate temperature. Indeed, FIDAP can demonstrate 
recirculations a) as the jet emerges from the nozzle due to flow 
separation, and b) on the sides of the hot susceptor due to bouyan- 
cy-induced convection. These recirculations are not detrimental 
because from the photograph one can see that particles trapped in 
the recirculation regions do not diffuse into the reagent jet, and 
we estimate that reaction product vapor species are also unable to 
penetrate appreciably into the jet. The Richardson number based on 
nozzle to substrate distance is -0.07 for this case. Lower flow 
rate (Re<100) measurements at higher substrate temperatures are 
hampered by buoyancy effects in the jet above the substrate. We 
discuss elsewhere [5] the operating conditions needed to avoid 
bouyancy— induced convection in such reactors. Our numerical 
analysis is capable of handling such non-ideal behavior. 

For the results discussed below we fix the gas flow rate at 5160 
seem (Re=500 for N 2 ) and the substrate-temperature-based TTIP 
concentration at C =2. 5x10^ Kgmol/m 3 . Deposition rates are reported 
as effective reaction probabilities e defined as e=n" / [ (1/4) vC s ] , 
where n" and v are the molar flux and mean thermal speed of TTIP, 
respectively, evaluated at the surface. 

Figure 4 depicts Soret diffusion and temperature-dependent gas 
transport— property effects on predicted rates at the stagnation 
point for transport controlled conditions. The rate reduction due 
to Soret transport is >75%, whereas the effect of variable 
properties is <15%, for T>1000K. 
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The model predicts reasonable deposition rates over the entire 
regime from surface kinetics to transport control by incorporating 
the surface kinetics extracted from our experiments (Fig. 5). With 
the gas-phase kinetics of [4], the model agreement with experiment 
is less satisfactory, though with similar trends. As the depletion 
of TTIP by the gas-phase reaction increases (i.e. chemical sublayer 
gets thicker) at higher substrate temperatures, thereby decreasing 
the TTIP surface flux (CVD rate) , FIDAP 's handling of Soret trans- 
port becomes more inconsistent. This is because the constant D T , 
used by FIDAP as defined above for the calculation of Soret trans- 
port, does not exhibit the expected diminishing Soret contribution 
as TTIP concentration diminishes. Hence, our rate predictions with 
Soret effect above substrate temperatures of 1200K are not shown. 
The cause (s) of our underprediction of the steepness of the high 
temperature rate fall-off is currently under investigation. 

The onset of homogeneous nucleation is inferred from a sudden 
drop in deposition rate at a certain substrate temperature. The 
reliability of the asymptotic theory prediction of this temperature 
depends on whether the chemical sublayer is indeed "much" thinner 
than the thermal boundary layer at the prevailing conditions. 
Figure 6 depicts the TTIP mass fraction profiles calculated by 
FIDAP at the stagnation point at a substrate temperature of 1200K 
with and without gas-phase reaction. We also plot the difference 
between the two mass fraction curves to demonstrate the extent of 
gas-phase reaction, as well as the corresponding temperature pro- 
file. The predicted rate at this temperature with gas-phase che- 
mistry is observably lower than the one without gas— phase chemist- 
ry. If the peak of this curve is taken as a measure of the thick- 
ness of the chemical sublayer, it is substantial (30-s) compared 
with the thermal boundary layer thickness. The amplitude and width 
of the difference curve would be larger at even higher substrate 
temperatures. Therefore, under such conditions predictions of the 
full numerical approach are expected to be more reliable than 
predictions of the simple asymptotic theory. Also noteworthy is 
the slight "enrichment" in the mass fraction profiles away from the 
surface due to the Soret effect. 

Figure 7 shows that the predicted rates on the susceptor are not 
uniform in the radial direction beyond the substrate radius (6.35mm 
vs. 16mm for the susceptor) . However, because the deposition rate 
is roughly uniform across the substrate, we have compared rate 
measurements to stagnation point rate predictions (Figs. 4 & 5) . 

5. CONCLUSIONS 

Our numerical model can correctly describe the flow field in our 
experimental atmospheric pressure, impinging jet reactor used for 
studying the CVD and gas-phase nucleation of titania. The model 
points to some buoyancy and flow recirculation effects inside the 
hot substrate reactor affecting the "ideal" stagnation point flow 
behavior. It can guide the interpretation of deposition rate mea- 
surements over a substrate temperature range (600-1400K) covering 
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surface-kinetics-, gas-phase transport-, and gas-phase reaction- 
governed regimes. Soret transport effects on predicted rates [6] 
are significant (-75% reduction at 1000K) , but variable property- 
effects are modest (<15% ) . Our numerical model extends the capabi- 
lity of the earlier asymptotic theory to predict the onset temper- 
ature of nucleation and the associated reduced deposition rates to 
include conditions where the chemical sublayer is not "much" thin- 
ner than the thermal boundary layer. Future work will include 
effects of flow rate and TTIP concentration [7], and a more careful 
examination of the high temperature rate fall-off regime. 
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To exhaust and scrubber 

Figure 1. — Schematic of the experimental reactor, after [1]. 
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Figure 2.— The computational domain. All dimensions in centimeters, cm 
and angles in degrees. 
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Figure 3. — Flow field for argon at Re = 750; T substrate = 900 K - ( a ) FIDAP simulation (gas velocity vector/ 
streamline contours), (b) Experimental visualization with titania seed particles. 
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Figure 4. — Predicted deposition rates with diffusion limitations showing variable property and Soret effects. Experimental 
data shown for reference. 
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Figure 5. — Comparison of predicted deposition rates (including both gas and surface kinetics) with experimental data. 
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Distance from substrate, cm 


Figure 6. — The effect of gas phase reaction on mass fraction profile at 1200 K. Soret effect is included. 
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Figure 7. — Radial variation of predicted deposition rates across susceptor surface. 
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